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ANDREI A. KAPAEV 

Abstract. For the bi-orthogonal polynomials with the third degree polyno- 
mial potential functions, the 3x3 matrix Riemann-Hilbert problem is explicitly 
constructed. The developed approach admits an extension to the bi-orthogonal 
polynomials with arbitrary polynomial potentials. 



1. Introduction 

The classical asymptotic theory of the orthogonal polynomials |) extended 
to the polynomials orthogonal w.r.t. the weight e~ NV ^ || [|, [|, || has gained an 
essential progress after introduction of the Riemann-Hilbert (RH) problem approach 
and the steepest descent method ||. A particular implementation of these 
methods to the orthogonal polynomials on the real line can by found in |5| a 
similar study of the orthogonal polynomials on the circle is given in . 

Further extensions of the notion of the orthogonal polynomials motivated by a 
number of applications to a random matrix theory, integrable systems, approxima- 
tion theory and combinatorics include the generalized orthogonal polynomials and 
the bi-orthogonal polynomials. In the former case, the sequence of the polynomials 
is orthogonal w.r.t. a sequence of measures Jl2| , 

J P n (X)dp m (X) = S nm , 

while in the latter case, two sequences of the polynomials are orthogonal to each 
other w.r.t. a two-dimensional measure, see [[u], [TBI pi , 

[ [ Pn(X)Q m (OM^O=Snm- 

Jr Jr 

In the simplest case, the two-dimensional measure has the form of the product 
d(i(\,£) = e -v(A)-w(0+A5 dXd £ where the p i ynom i a i s y(X) and W(0 are called 

potentials. Integration over £ in the latter two-fold integral yields a sequence of 
measures, dp m (X) — <f) m (X)dX = J R Q m (£,) d[i(X, £). The functions <p m (X) here are 
called the dual functions, see |l5j . 

The generalized and bi-orthogonal polynomials are associated to a completely 
integrable system coming from the ^-deformations and the Virasoro constraints and 
describing typically the reductions of the 2-Toda lattice, sec |l2|, |l3|, [IB], M. 
However, in spite of the algebraic properties of the generalized and bi-orthogonal 
polynomials are well known, the knowledge of the asymptotic properties of such 
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polynomials is very limited. Basically, this is due to the absence of an adequate 
formulation of the relevant RH problem. Indeed, in spite of the 2x2 matrix RH 
problem formulation for the conventional orthogonal polynomials Q admits a direct 
2x2 extension for the generalized orthogonal polynomials [[l2| , the relevant version 
for the bi-orthogonal polynomials E3 shows its non-local nature. 



Fortunately, the study |15 reveals the isomonodromy structure associated to the 
bi-orthogonal polynomials and hence the principal possibility to formulate annxn 
matrix RH problem which would enjoy the properties similar to those for the 2x2 
matrix RH problem for the conventional orthogonal polynomials. In what follows, 
we construct the 3x3 matrix RH problem for the bi-orthogonal polynomials with 
the third degree polynomial potentials, as well as the 3x3 RH problem for the 
relevant dual functions. In spite of its less physical importance, this case provides 
us with the opportunity to develop the technique in the simplest non-trivial case. 
Indeed, the bi-orthogonal polynomials for the second degree potentials are reduced 
to the classical Hermite polynomials jl3]. In some more involved case deg V(X) > 2, 
deg W(£) = 2, the bi-orthogonal polynomials can be expressed in terms of the semi- 
classical orthogonal polynomials related to the 2x2 matrix RH problem studied in 
@ |, and other papers. 

We stress that the method presented below can be extended to arbitrary poly- 
nomial potentials. We also point out a particular importance of the paper [ fl7| 
which is useful for construction and justification of the RH problem for the dual 
functions for the class of the weights with the rational log-derivatives. After the 
original version of the present paper was posted to the internet, the Montreal group 
presented their methodology for construction similar RH problem ]l8|] . Their idea 
for evaluation of the RH problem for the dual functions based on the study of the 
path integrals coincides with the presented below but, in contrast to our cubic case, 
they consider the arbitrary polynomial potentials. As to the RH problem for the 
original (wave) function, the authors of Q rely on the so-called "duality pairing" , 
while our approach is based on the explicit integral representation for the wave 
function. 

The paper is organized as follows. In Section]^, we recall the matrix differential 
and difference equations satisfied by the bi-orthogonal polynomials and their dual 
functions. In Section ^, we construct particular matrix solutions for the differential- 
difference syste ms for t h e dual functi ons whi ch gives rise to a RH problems for the 



dual functions ( 3.3.11 )-( 3.3.13|) and ( 3.3.15 )— ( 3.3.17 ). In Section || we construct 



particular solutions of the differential-dif ference s ystems for th e bi-or th ogonal poly- 
nomials and the relevant RH problems (|4.4.32|) - (^.4.34D and (|4.4.35|) - (|l.4.37| ). 



2. Equations for the bi-orthogonal polynomials 

Below, we consider the monic polynomials p„(A), g m (£) satisfying the orthogo- 
nality condition 

' d\ [ dt; Pn (\)q m (t;)e- VW - w< t )+tx Z =hl5 n>m , (2.2.1) 

71 ^72 

where 

^(A) = iA 3 +xA, W{0 = \e + Vt, (2.2.2) 
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x, y, t G C, t 7^ 0; the contours ji, i = 1,2, are the complex linear combinations of 
the elementary contours, 

2 

fff eC, 1 = 1,2, (2.2.3) 

3=0 

where each Tj is the sum of two rays, 

r^- = (e^^-^oo.O] U [0,e iS rU-Voo), j = 0,1,2. (2.2.4) 

Because T + Ti + T 2 = 0, one of the parameters g^ (resp., g^) can be put to 

zero; one of the nontrivial parameters g^ can be normalized to the unit. Thus the 
set of contours jj and therefore the set of the monic bi-orthogonal polynomials is 
parameterized by three constant complex parameters^} 
Introduce the wave functions 

^n(A) = ^p n (X)e- v ^, «MO = -L qm (£)e- W ^, (2.2.5) 
tin f^tn 

and define the dual functions ^Vi(£), 4>m{X) which are the Fourier-Laplace images 
of the wave functions p| , 

Mt)= f ^n(A)e* A «rfA, 4>m(X)= [ cf> m (0e m <%- (2-2.6) 



The orthogonality condition (2.2.1) now reads 



dX ^^(A)0 m (C)e tA5 = / V«(A>MA)dA= / M04>m{t) = S v 

7i J 72 J 7i *^ 72 



(2.2.7) 



It implies certain relations between the introduced functions ( 2.2.5|) , (2.2.6). We 



refer |l^j for the general case and present here the final result for our particular 
situation: 

n+l n+2 

AV>„(A) = ^2 a n ,m^m{X), d\ip n (\) = -i b, hm ip m (X), 

m=n—2 m—n—1 
n+l n 
d x ip n (X) = ^ u n,mi>m{X), d y i/j„(X) = - ^ V n , m \j) m (A) , (2.2.8) 
m— n 

n-1 

d t 1pn(X) = W n 1pn(X) - ^ A n,m1pm{X), 



m—n—1 



m— n — 3 



1 More general parameterization of the set of contours is introduced in ]is| , however, in our 
cubic case, both kinds of parameterizations are equivalent 
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m+1 



m+2 



n=m— 2 n— m— 1 

m+1 m 

dy<t>m(£,) = ^2 V n,™<t>n{£)-> Qe^m(f) = ~ ^ ^n,m</>n(£)> (2.2.9) 
n—m n—rn—1 
m — 1 



n—m— 3 



where 



&n,n+l — ^n+l,n 








h n +i 




, h n +2 

On.n+2 — , , 

tn„ 


— "n,n+l — 


— V n +l,n — 


h n 


Q-n+2,n — 




In /i n , 


^n,n 


= d y ln/i n , 


A 


—2w n = 


2d t ln h n ; 






' Y^m+1 


_2 ^n,k^k : Tn ? 


< n 


- m < 3, 




771 * 


y-vn+1 

n 2-^k—m 


_2 ^TL^k^k^m ■> 


-3 < n 


- TO < 





Using the definition fl2.2.6|) and the equations above, it is straightforward that the 
dual functions satisfy 

n+l n+2 



m=n—2 m—n—1 

n+1 n 

m—n m=n—l 
ra+3 

m=n-\-l 

m+1 m+2 

d\<fr m (\) = tb njm 4> n (\), A$ m (A) = a n , m 4> n (\), 



(2.2.10) 



n—m— 1 



(2.2.11) 



n—m — 2 

m+1 m 
d y <j) m {\) = ^ ^n,m0n(A), d x (j) m {\) = - ^ ^n,m0n(A), 
n—m n=m— 1 

m+3 

dt<f>m(\) = -u> m 0m(A) + ^4 n5m (^ n (A). 

n=m+l 

The above expressions can be written in the matrix form [ Ejj . Indeed, the 3- 

vector \P„(A) = (■?/!„ (A), V'n-i(A), V'n-aM) satisfies the system of difference and 
differential equations, 



*„+i(A) = i?„(A)*„(A), ^(A) = A„(A)f„(A), 
<M "" =C/ n (A)*„, ^ = K(A)* n , ^ = ^„(A)* ? 



(2.2.12) 



d x 



dy 



where 



' X—Q, n ,n a nn — X 0,n,n — 2 s 

On, n+l a-n,n + l Q>n,n+1 

1 
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A n (X) = ^-Dn,n+2-Rn+l(A)-R„(A) — D„, n+ ii?„(A) — — D n ,n- 1 1 (A) , 

Ai,m = * diag(6„ !m , 6 n _i jTO _i, & n -2,ro-2) , 
f/ n (A) = £/n,n+l-Rn(A) + ?7 n ,nj f^n.m = diag(w n , TO , ti n -l,m-l) U n - 2 ,m-2) , 
V n (X) = —V nin — V n , n -\R n \(\), V n ,m = diag(«„ iTO , V n -i tTn -i , V n -2,m-2) , 

W n {\) = W n , n - An^R-^iX) - A n , n - 2 R-1 2 (X)R-^W - 

W„, n = diag(w„, w„_i, w n - 2 ), A n , m = diag(A„, m , A n _i >m _i, A„_ 2 , m - 2 ) ■ 
Similarly, the 3-vector <J> m (£) = (0 m (£)> m _i(£), 4>m- 2 {€)) T satisfies the system 

*m+l(0 = Q m (0* m (0. ^(0 = Bm(0* m (0. 

9$ 9$ 9$ (2.2.13) 

ox ay at 

where 

f£-b m , m b m -!, m 6 m _ 2 , 



f>m + l, m f>m + l,m &m+l, 

Qm(0 = | 1 

o i o 

B m (0 = —D m + 2 ,mQm+l{€)Qm{€) — ^m+1 - Vl,m<3 m -l(6i 

T^m,n = t diag(a TOi „, a TO _i j7l _i, a m - 2 ,n- 2 ) , 

W m (£) = ~U rn , m — U m — l,roQ TO _ Vm(£) = Ki+l,mQm(() + ^m,mj 

-^m-3, m Q- 1 _3(0^ 1 -2(0Q™-l(0- 

The matrix equations for the 3-vector of dual functions, 

*n(0 = (^n(0,^n-l(0,^n-2(0) T , 

are as follows, 

^(0 = ^n(0*n(0. ^(0 = K(0*n(0. ^(0 = W" n (0*n(0> 

(2.2.14) 

where 

&n-l,n £ — &n-l,n-l 



&n-l,n + l &n — l,n + l &n — l,n + l 

#n(£) = ( 1 

1 

in(0 = ©n,n+l^(0 + ^n,n +X>„,„_l^i: 1 (0 +©„,n-2^ 2 (0^-l(0. 
Un(0 = U n , n+ lRn(0 + Un.n, V n (0 = ~V n , n - V n>n -! R^i 1 (0 , 

W"n(0 = -W n , n + A n , n+3 R n+2 (0R n+1 (0R n (® + 

+ A n ,n+ 2 Rn+l(£,)Rn{£,) + A n ,n+lRn(0- 
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In the very same way, the 3- vector 
satisfies the system of equations 



$m+l(A) = Q m (A)$ m (A), 



■(A) = B m (X)$ m (X), 



! " Kl (X) =Z4,(A)* m (A), ^(A)=V m (A)$ ro (A), ^(A) = W m (A)$ m (A), 



dx 
where 



dy 



dt 



(2.2.15) 



Q-m.,m-l A — a m „i iT7T _i a7Tl-2,rJ^-l , 



»(A) 



Om + l,m-l Clra+l.ra — 1 Om+l,m — 1 

1 



1 

B m (X) — Dm+l,mQm(X) + D mm + £*m-l,mQ m _i(A) + D m —2,mQ m —2 

(A)Q m 1 _ 1 (A), 

Um(X) = —Um,m ~ U m —l jm Q m _ 1 (X), V m (A) = V m +\ (A) + y m 



Vv m (A) 



(A)Q m+ i(A)Q m (A) 



.4 



(A)Q m (A) + ^4 m +l,mQm(A). 



The compatibility conditions of the above equations yield a reduction of the 
2-Toda lattice, see |12}. On the other hand, this nonlinear system describes the 
isomonodromy deformations of the 3x3 m a trix diff e rential equat ions in A and £. 
The matrix solutions of the systems (|2.2.12| ), (|2.2.13|) , fl2.2.14|) and (|2.2.15| ) give rise 
to a matrix RH problems whose data can be found via the asymptotic analysis of the 
A- and ^-equations with the polynomial matrices A n (X), B m (£) and A n (X), B m (£) 
using the WKB technique, see |1S| ]. It is crucial however, that it is possible to avoid 
the WKB analysis expanding the basic idea of 0. In what follows, we construct 
certain particular solutions of the above systems in terms of the (unknown) bi- 
orthogo n al polyn o mials u si ng mere ly the fact of the existence of the linear systems 
( |2.2.12|) , (12.2.131) , Q2.2.14D , ( |2.2.15|) rather than the systems themselves. 



3. The particular solutions of the matrix equations and the 
rlemann-hllbert problem for the dual functions 

Because of the mentioned above independence of the relevant monodromy data 
on the deformation parameter t, without loss of generality, we restrict ourselves to 
t > 0. This assumption allows us to simplify our calculations while the final result 
will be valid for arbitrary t E C\{0}. 

Let us introduce the auxiliary functions 

$°(0=/ ^„(A)e tA «rfA, $<£{X)= ( m (Oe tA «^, 



(3.3.1) 

where the contours 1^, j = 0, 1,2, are defined in @.2.4| ). Due to ( |2.2.6| ) and ( |2.2.3[ ), 
the dual and auxiliary functions are related to each other via 

MO = X^W 3 ®, MA) = j^gf } ^(A), 
j=o j=o 



(3.3.2) 
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and satisfy the same equations (2.2.1C), (2.2.11). Taking into account the orthog- 
onality condition (2.2.7), we find that the functions 



2m 



2m 



C-A 



dC, 

(3.3.3) 



for n, m > 2 also satisfy the same equations. 

Observing that $n , j = 1, 2, (respectively, (j)m\ j = 1,2) ( |3 .3.1 ) are combina- 
tions of the independent Airy functions and their derivatives, 



dX 



dm (A) 



-q m (d T ) 



e 3 



dl; 



T=t£ — X 

r=tA-s/' 



(3.3.4) 



we construct the particular piece- wise holomorphic solutions of the system of the 
3x3 matrix equations (2.2.14) and (2.2.15), 



*n(£) 



and 



$m(A) 



F n _^) 
WnUO ^nUi) Fn-2(0y 



/ ^ } (A) 0m' (A) G m (A) 

0m-l(A) ^Li(A) G m -l(A) 

^-2 (A) G m _ 2 (A), 



n > 4, 



m > 4. 



The jump property of the Cauchy integral yields the relations 



^ (0 - F~(0 = {gf - g^JMt), arg^ = - 1), j = 0, 1, 2 



2?r 



G m (A)-G m (A) = (.g J (1) - 5 W 1 )0 m (A), argA=^(i-l), j = 0,1,2, 



2tt 



(3.3.5) 



(3.3.6) 



(3.3.7) 



where = g$ . Thus we find that the matrix functions ^n(£) and <& m (A) have 
the following jumps across the rays £j ={(eC: & r g£ = x0 — ■*■)}' ' = ^' ^' 
oriented towards infinity, 



*£(0 = *»tt)|o i (4 2) -^i)^-^) 

^0 1 

(l {g? - gflM ] - a^Y 
**(A) = *«(A) 1 {g? - gfl){$> - 
\0 1 



£e£„ j = 0,i,2, 



Ae4 J = 0,l,2. 



(3.3. 
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Using the well known asymptotics of the Airy integrals in the complex domain, 



;A 3 +tA _ 



-iV^r-^e-t^^l + OCr- 3 / 2 )), argre, 



dA = ? :^r- 1 / 4 e-^ 7 (1 + G(r- 3 / 2 )), argrG (f ,2ir), 



e 3 



A 3 +tA 



dX 



r 2 



(3.3.9) 



(3.3.10) 



the asymptotics of the Cauchy integrals for ( 3.3. 3| ) , 

F n (0 = -^- i C n - 1 e W ^(l + O(C 1 )), 
G m (A) = -^A- m - 1 e^ A )(l + 0(A- 1 )), 

we construct the RH problem for the dual functions ip n {£,)- 

Riemann-Hilbert problem 1. Find a piece-wise holomorphic 3x3 matrix func- 
tion \&„ (£) wifft i/ie following properties: 
1. As £ — + oo, 



/ef(*?) 3/2 — (*« 1/2 



2-n-i S 



27TI 



V 



eif 3 +f« 



(3.3.11) 

i?. Across the rays arg£ = ^-(j — 1), j = 1,2,3, oriented towards infinity, 



^n H (£) h- as the jumps 



,RH 



2tt 



(3.3.12) 



where plus and minus indicate the limiting values of '(£,) on the jump contour 
from the left and from the right, respectively, and 



Si = 



1 (g? - $)<gV - gP) 

-i i ii^-a^-^l' 



o o 



S 2 = 



5 3 = 



i -i (g? -gifW?-^) 

,0 1 

1 ( 9 P - 9?%^ - 9 P) 

-i 1 i(£> - gPXgP - tf) 







1 
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and across the ray arg£ 
tion holds, 



■| oriented towards infinity, the following jump condi- 



RH-\ 



arg£=--, S = 




(3.3.13) 



The dual functions 4> n (£), ip n -i(£) and tp n -2(0 form the vector 9 n (Q related 
to the solution of the RH problem [j] by the following equations, 

(*™(ORW, arg£e(-§,0), 



*n(0 



[*™(OS^S^S^RU, 



arg£e (0,f ), 
arg£e(f,f), 



(3.3.14) 



where 



V 




Riemann-Hilbert problem 2. find a piece-wise holomorphic 3x3 matrix func- 
tion ^>^ H (X) with the following properties: 



1. As A 



00, 



.(tA)^-i 7 

(tA^-j^A) 1 / 2 






jg(-ir(*A)9-i 




/t m \ — m— l 

27T'i 


£-(-!)— i(tA) 1 ^ 


1 
4 






1 
4 


hm-2 \-m+l 

27TI 













), 





1 

C 3 


A 3 +:rA / 



(3.3.15) 



2. Across the rays argA = — 1), j = 0,1,2, oriented towards infinity, 



$m H {\) has the j 



jumps 



2tt 



<f£ K+ (A) = argA=-0'-l) 



(3.3.16) 



where plus and minus indicate the left and right limits of $^^(A) on the jump 
contour, and 



Ti = 



Is = 












1 






-^ 2) ) 







i 




— i 








1 






-^ 2) ) 







i 







(^ 




-«?>) 


1 













i 
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and across the ray arg A 
tion holds, 



TT oriented towards infinity, the following jump condi- 



3> 



RH- 



(A) = $^-(A)£, arg A 



TT 

3' 




(3.3.17) 



The dual functions m (A), m _i(A) and m _2(A) form the vector $ m (A) related 
to the solution of the RH problem @ by the equations 



where 



* m (A) 



>r(A)i? (2) , 

$r(A)Tr^ (2) , 

^(A^ffW, 
&* H (\)f 3 -%-%-iRW 



argAG(-f,0), 
argAe (0,f ), 
argAG (f ,f ), 
argAe (f ,f ), 



/ (2) 
' 01 ■ 

.(2) 



(2) 
,(2)> 







(3.3.18) 



4. The particular solutions of the matrix equations and the 
rlemann-hllbert problem for the wave functions 

The method of construction of the Riemann-Hilbert problem for the wave func- 
tions ^n(A) and <fim(£) is more involved because of less elementary structure of their 
integral representations in compare to those for dual functions. 

Let be an oriented contour connecting a finite point £o (or Ao) with infinity 
within the sector of the exponential decay of the function ipn (£) ( |3.3.4 ) (resp., of 
</>m^(A)). Namely, let £q^ be asymptotic to the ray 



[0, e 3 r - J oo). 



Let Tj be an infinite oriented contour asymptotic to the rays 



(e 



2 -rU-i) r 



oo,0] U [0,e^ (j '"5) O o) J j = 0,1,2, 



located within the sector (4.4.2) and intersecting the ray 



ij = [0, 



oo), 



(4.4.1) 
(4.4.2) 
(4.4.3) 



so that Tj n ij = {£,-} (resp., Tj D ij = {Xj}). 

The "bricks" which we use to build up the appropriate integral representations 
are the "inverse" Laplace-Fourier transforms 

t 



t 



Til lii) 



and 



(A)dA, G m (0 



7TJ Jf 



F„(0e" tA «^, (4.4.4) 



G m (A)e ^ cZA, (4.4.5) 



The functions (4.4.4) satisfy the differential equations in A, x and y in (2.2. 



but not the recursion relation and the differential equation in t because of the 
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appearance of the inappropriate off- integral terms after integration by parts. The 
similar observation holds true for (4.4.5). More in details, 



A^ fc) (A) = — e- tA «°^ fe) (£o) + appropriate 



terms, 



(4.4.6) 



A^n(A) = ~ e ~ tXij ( F n(£j) ~ F n(€jj) + appropriate terms = 

1 2 
= -e~ tx ^ {gf - g\%) ]T g^ffi fo) + appropriate terms. 



(4.4.7) 



k=0 



In second line of (4.4.7), we have used the jump condition (3.3.7) and the definition 
(3.3.2). Combining (4.4.C) at £o = £ 3 and (4.4.7), it is possible to eliminate the 
off- integral terms and find such a combination F$'(X) of ^^(A) and F n (\) which 
satisfies the system (2.2.8), 



(4.4.8) 



k=0 



Using identity Ylj^HO = 0, it is possible to eliminate one of the values ^i fc ^(^) 
from ( 4.4.7 ) and therefore the respective term in ( 4.4. 8 ). By technical reas ons, w e 
prefer to use the latter possibility and introduce the following solutions of ( 2.2.8| ): 

^ 0) (A) = F n (X) - ( 9 P - g? - 9^W{X) - (9^ - tfW^ - gPtfPW, 

(4.4.9) 

WW = F n (X) - (g[V £ ))(,« _ 3 W)^)(A) - fo?> - gf){g^ - ^)^(A), 

(4.4.40) 

^ 2) (A) = F n (X) - (£> - g^ - g^HX) - {gf - g^)(g^ - ^)^(A), 

(4.4.44) 

In the very same way, 

(4.4.42) 

G£>(0 = G m (t) - (ffW - - ^ 2) )0W(e) - - - g^W(0, 

(4.4.43) 



GW(0=G m (0-(fl 



(4.4.44) 

Define the 3x3 matrix functions 

/ Vn(A) i^ 0) (A) Fi 1} (A) 
V„_i(A) ff\(A) C\(A) 
VV»- 2 (A) Ff 2 (A) F« 3 (A) y 



* n (A) 



n > 4, 



(4.4.45) 



m>4. (4.4.16) 
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and 

/ MO G$(0 G#(0 \ 
*m(0 = Ml(0 Gf^C) ^^(0 

Vm 2 (o g^_ 2 (o G« 2 (oy 

The asymptotics at infinity of ?/> rl (A) and </>,„(£) is elementary, 

IMA) - ^-e^ A ' - A (l + ©(A- 1 )), MO = |_ e -i« -»«(! + OiC 1 ))- 
n n ii m 

(4.4.17) 

The asymptotics of Pi j) (A) and Gm (0 can be found using the conventional steepest 
descent method, 

F(D(A) = ^L(tA)- 2 * 1 -ie-i( tA ) s/2 +''( tA ) 1/a (l + 0(A- 1 /2)) + 

+ 2(ff« - ^ 2) )( 5 « - ^))^e-i A3 - A (l + 0(A -1 )), arg A e (-f , 0), 

(4.4.18) 

F«(A) = |^(tA)- 2 * i -ie-i(* A ) 8/a +»(* A ) 1/a (l + 0(A- 1 /2)) + 
+ 2(. 9 f ) - ff f)(. 9 «- 5 «)^ e -^ 3 -^(l + 0(A- 1 )), arg AG (0, f ), (4.4.19) 

+ 2(g^-g?%i 1) ~g[ 1) )^e-i x3 -* x (l + <D(\- 1 )), argAe (0, f ), (4.4.20) 
Ff(A) = -^(-l)"(a)- £ * 1 -Jei(* A ) 3/2 -^* A ) 1/2 (l + 0(A- 1 /2)) + 



2tt 3 / 2 

\ n 

'X^ - S^jJ-e-^'-^l + ©(A- 1 )), arg A G (f , f ), 

(4.4.21) 



and 



Fp)(A) = -|^(a)-^-i e -i( tA ) 3/2 +^ A ) 1/2 (l + G(A- 1 / 2 )) + 



+ 2( ff f - .g^)(.g« - ^>)£e-** a -*(l + ©(A -1 )), arg A e (f , f ), 



(4.4.22) 



^ 2) W = -^(tA)-^-i e -t^) 3/2 +^) 1/2 (1 + 0(A -1 / 2 )) + 



2tt 3 / 2 

W - ^^e"^ 3 -^^ + ©(A- 1 )), arg A e (f ,2tt). 

(4.4.23) 
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For the asymptotics of G m (Q, we have 

<^ } (0 = ^(io-^-e-i^ 3/2+ ^) 1/2 (i+o(r 1/2 )) + 



+ 2( ff « -&>)(g[*> -^jfe-^-^l + Or 1 )), arg^e (-f ,0) 

(4.4.24) 



^(0 = ^(i0-^- e -i^ 3/2 +^) 1/2 (i+o(r 1/2 )) + 

+ 2( 5 «-^ 1) )( 5 f-^ 2) )|^ e -^ 3 -^(l + 0(r 1 )), arg^G (0, f ) (4.4.25) 



^ } (0 = -|^(-in^r^- ^ (t5)3/2 -^ )1/2 (i + o(r 1/2 )) + 

+ 2( 5 « - g^){g^ g^)i-e-^-yi(l + ©(r 1 )), arg£ £ (0, f ), (4.4.26) 

^(0 = -^(-i)™(io-^-ei^ 3/2 -^) 1/2 (i + o(r 1/2 )) + 

+ 2(. 9 « - gV)^ g ?)^e-^-yt(l + OiC 1 )), argC € (f , f ), 

(4.4.27) 

and 

+ 2(.4 1} - ^(sf - gP)i-e-&-«(l + OiC 1 )), argC € (f , f ), 

(4.4.28) 

^(fl - -|^(i0-^-e-i(«) 3/2 +^*« 1/2 (i + o(c 1/2 )) + 

+ 2( 5 « - gV)(g[V g ^ e -H 3 -v^i + QiC 1 )), arg£ e (f ,2tt). 

(4.4.29) 

Using the above asymptotics and the linear constraints for F„ \ Gm\ 
^ 0) (A) + PP(X) + F^(\) = 2 ffF ^„(A), 

9F = - 9[ 1] ) + ffPW - + £W 9 ( 2 \ (4.4.30) 

and 

g^(0 + gW(0 + (5^(0 = 2^(0, 

9G = W ~ 9?) + W ~ ^) + flW - <?f), (4.4.31) 
we find the RH problems for our bi-orthogonal polynomials. 
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Riemann-Hilbert problem 3. Find a piece-wise holomorphic 3x3 matrix func- 
tion ^^ H (X) with the following properties: 
1. As A — ► oo, 



*r(A) 



( £ -^(-i)"(u)-^- 3 ^ ( u)-^- 3 



e -|A 3 -xA o \ 

e §(tA)V»-,(tA)^ o 

e -i(t\f'*+y(t\)^ j 



(4.4.32) 



2. Across the rays argA = -^(j — 1), j = 1,2,3, oriented towards infinity, 
^n H W has the jumps 

*2 H+ W = *™-(\)S j , argA= y (i-l), (4.4.33) 

where plus and minus indicate the limiting values of ^^ H (\) on the jump contour 
from the left and from the right, respectively, and 

A - tfW - tf) 2{g? - £>){gV - gP? 

#1= 1 

\0 -1 1 



S 2 = | 1 1 

v 1 



1 

v -1 1 

and across the ray arg A = — J oriented towards infinity, the following jump condi- 
tion holds, 



*^+(A) = *^-(A)S, argA = -|, E = | (4.4.34) 



'10 
0-1 
,0 1 



The wave functions i/) n (\), Vn-i(A) and ^>„_2(A) are just entries of the first 
column of $>% H (\). 

Riemann-Hilbert problem 4. Find a piece-wise holomorphic 3x3 matrix func- 
tion Q 1 ^ 1 (£) wit/i i/ie following properties: 
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1. As £ — ► 00, 



th„ 



+ i i 

2 4 



ll 



-2^(-i) m (*e) 
.^ ( _ ir -i (t0 ---i 

|^(-i r -2 ( ^)-^ 



2 



o 

V 



es 1 



I_I ith m -2 (j.C\- 

4 "2^372- 




e -l(te) 3/2 +-w) 1/2 



n - 1 1 

2 4 



(4.4.35) 



2. Across i/ie rays arg£ = ^(j — 1), j = 1,2,3, oriented towards infinity, 
®m H (0 has the jumps 



(4.4.36) 



w/iere pfozs and minus indicate the limiting values of <&^ H (£) on the jump contour 
from the left and from the right, respectively, and 



Ti = [ 



-1 



1 
vO 

(1) (2) J2) 



T 3 = | 
^0 



1 
1 

(1) Jl)w J2) (2)s 



1 

-1 



and across the ray arg A = — ^ oriented towards infinity, the following jump condi- 
tion holds, 



$r + (C) = $r-(OS, arg£ 



3' 



'10 
v = | -1 I 

,0 1 



(4.4.37) 



The wave functions <j> n (\), n _i(A) and n _2(A) arc just entries of the first 
column of $% H (t;). 
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